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Event Shapes and Power Corrections in ep DIS 
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Deep-inelastic ep scattering data, taken with the HI detector at HERA, are used to study the event shape 
variables thrust, jet broadening, jet mass, C parameter, and two kinds of differential two-jet rates over a large 
range of "relevant energy" Q between 7GeV and 100 GeV. 

The Q dependence of the mean values is fit to second order calculations of perturbative QCD applying power 
law corrections proportional to 1/Q^ to account for hadronization effects. The concept of these power corrections 
is tested by a systematic investigation in terms of a non-perturbative parameter Qp_i and the strong coupling 
constant. 



1. INTRODUCTION 

Due to asymptotic freedom perturbative QCD 
is a powerful theoretical tool for the investiga- 
tion of strong interactions at high energies Q. 
Given, however, the precision of today's measure- 
ments over a large range of energy, deviations 
diminishing proportional to 1/Q^ are observed. 
Depending on the sensitivity of the observable 
on long-distance, non-perturbative effects, these 
power corrections may be sizable. Less inclu- 
sive quantities than e.g. the total cross section 
of e"'"e~ —^ hadrons, where p = 4, exhibit much 
larger corrections; in case of event shapes they 
are typically proportional to l/Q- 

One possibility to account for the non- 
perturbative formation of hadrons out of quarks 
and giuons is to invoke phenomenological models 
implemented in Monte Carlo generators. As a 
drawback, one has to deal with several models, 
each again comprising a large number of tunable 
parameters. Alternatively, analytical formulae ei- 
ther inspired by a simple tube model [ |l[] or by a 
recent approach initiated by Dokshitzer and Web- 
ber [ m may be employed. A first study of event 
shapes and analytically parameterized power cor- 
rections in ep deep- inelastic scattering (DIS) was 
published in [ ^. The results presented here 
are based on the substantially improved and ex- 
tended analysis in [ H| . 



2. EVENT SHAPE VARIABLES 

The event shapes investigated are two kinds of 
1-thrust, T and tc, the jet broadening B, the 
jet mass p, the C parameter, and two kinds of 
differential two-jet rates, y/j and ykt , with y 
denoting the transition value of an event from 
(2 -I- 1) — > (1 -t- 1) jets. All of these infrared and 
collinear safe variables are defined in the Breit 
frame of reference. Except for the two jet ob- 
servables y/j and j/^^, they exploit the property 
of this frame to maximally separate the current 
quark jet from the proton remnant by evaluating 
the current hemisphere alone. They are normal- 
ized to the total momentum resp. the total energy 
in the current hemisphere; the precise definitions 
may be found in [ H, H. In case of the y variables, 
the employed jet algorithms, i.e. the factorizable 
JADE and the kt algorithm [ H , implicitly isolate 
the proton remnant. 

The energy scale relevant in DIS is given by 
the four-momentum transfer q of the scattered 
lepton as Q := \/—q^ and ranges from 7 GeV 
up to 100 GeV in this study. Note that a 
second quantity, e.g. the Bj0rken scaling vari- 
able x:=—q^/2P-q with P being the four- 
momentum of the incoming proton, is necessary 
to fix the DIS kinematics. Here, x varies between 
w 10~3 and w 0.5 . 



3. POWER CORRECTIONS 

According to [ the energy dependence of any 
event shape mean (F) can be written as sum of a 
perturbative part (F)^'^'^ and a power suppressed 
contribution (F)^°^ . Identifying for simpUcity 
the renormaUzation scale /in with Q, the pertur- 
bative part, evaluated to O(a^) using DISENT [ 
is given by 

{Fr'' = c^Mx) c^siQ) + C2Ax) aliQ) . (1) 
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Note that in contrast to e"'"e~ annihilation parton 
density functions are necessary to describe the 
initial state in ep DIS. Apart from other quan- 
tities the coefficients therefore explicitly depend 
on X. 

Anticipating power corrections proportional to 
1/Q and 1/Q^ in the case of (j/fcj), the simple 
ansatz of (F)P°^ = X/Q resp. {F^"^ = n/Q^ is 
tested. However, fits of A or /i alone are not able 
to describe the data satisfactorily. Fitting also 
as{Mz) mostly leads to unacceptable values of 
as{Mz) owing to strong correlations. Triggered 
by the x dependence of (F)^*"^ , one could assume 
that (F)^°™ should be a function of x as well. 
Yet, for corrections with power p = I this is not 
expected [ p[ . 

This calls for the more elaborate approach of [ 
0], referring infrared divergences in the pertur- 
bative expansion to the power corrections. As- 
suming the existence of a well-behaved effective 
strong coupling even at scales below the infrared 
matching scale /i^ of w 2 GeV, one can deduce the 
relative size of the corrections V/Q to the means 
of different event shapes at the expense of a new 
non-perturbative parameter cEq: 



V = apM 



,16^ 



ao(M/) - as{Q) 



-^ fin ^ + 1.45 )a^(g) 



(2) 



where ap is a calculable F dependent coefficient. 
To resolve ambiguities a refinement to two-loop 
level leads to another global factor, the Milan fac- 
tor M, where TW' = 2A^/7r w 1.14 [|. 

As an example, fig. n] shows the unfolded dif- 
ferential distributions and a power correction fit 
for yfj, however, with ay^j = —1/4. 
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Figure 1. Left: Norm. diff. distr. of yfj cor- 
rected for detector effects. HI data (stat. un- 
cert. only) are compared with NLO calculations. 
Right: Corr. means of y/j (stat. and tot. uncert.) 
as a function of Q. The full line corresponds to a 
power correction fit with ay^j — —1/4. 

The results of the two-parameter fits to this ap- 
proach are presented in table ^ and fig. g. The 
universal non-perturbative parameter ao{fii = 
2 GeV) is confirmed to be « 0.5 ± 20% for the 
means of r, B, tc, p, and C, which are sizably af- 
fected by hadronization. The spread in as{Mz), 
however, is uncomfortably large. In contrast to 
earlier analyses of -B, a reasonable fit is obtained 
if the revised coefficient a'g [ [l^l is included. 

Both y observables exhibit small negative 
hadronization corrections. For yfj a power cor- 
rection with p ^ 1 and ay^j — 1 was proposed in 
the conference proceedings [ O], which has not 
been reconsidered in recent publications. From 
the emerging very low values for both, ao(/ij = 
2 GeV) and as{Mz), one can rule out ayj,j = 1. 
Turning the argument around, a reasonable fit 
can be achieved with ay.j = —1/4 while presup- 
posing the validity of eq. 0. For yk^ fits with p — 1 
do not work properly in accordance with the ex- 
pectation of a power p = 2. Unfortunately, the 
coefficients ap are not yet calculable for that case. 



4. SUMMARY AND CONCLUSIONS 

The improved and extended analysis of event 
shape means in ep DIS supports the concept of 
power corrections according to eq. 0. However, 
the uncomfortably large spread of values gained 
for as{Mz) leaves room for further progress. 



Table 1 

Results of two-parameter fits oiao{fij — 2 GeV) and as{Mz) according to the power correction approach 
initiated by Dokshitzer and Webber for the event shape means. The uncertainties are statistical and 
total systematic except for yfj where a coefficient Uy^j different from 1 is anticipated, k denotes the 
correlation coefficient between the two fit parameters. 
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PC 1/y^q7 + const, as suggested by Dokshitzer, Marchesini and Salam [ p_0[ 
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Figure 2. Results of two-parameter power correc- 
tion fits for the means of r, B, tc, p, and C in the 
form of la- and 2(T-contours in the (as, cfo)-plane 
including stat. and exp. syst. uncertainties. 

Like in e"'"e~ annihilation, it would be very in- 
teresting to compare to the power corrections ob- 
tained with event shape distributions. To that 
goal, resummed calculations are awaited for. 
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